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FREE  OSCILLATIONS  OF  AN  ATMOSPHERE  IN  WHICH  TEMPERATURE 
INCREASES  LINEARLY  WITH  HEIGHT 
By  C.  L.  Pekeris 


SUMMARY 


It  is  shown  that  when  the  temperature  in  the  atmosphere  increases 
linearly  with  height^  the  speed  of  propagation  of  long  waves  does  not 
approach  a  limit  with  increasing  wave  lengthy  as  in  the  case  of  an 
atmosphere  in  which  the  temperature  at  great  heights  is  assumed  to  be 
constant  or  decreasing,  but  increases  linearly  with  the  period.  The 
group  velocity  ultimately  also  increases  linearly  with  the  period  and 
becomes  equal  to  half  the  phase  velocity.  The  region  of  maximum  energy 
of  the  oscillation  is  shifted  to  increasingly  higher  elevations  as  the 
period  is  increased.  Whereas,  in  an  atmosphere  where  the  temperature 
gradient  at  great  heights  is  negative  or  zero,  the  tides  are  similar  to 
those  in  a  uniform  ocean  of  equivalent  depth  H  (about  8  km,  depending 
on  the  assumed  vertical  temperature  distribution),  the  superposition  of 
an  outer  envelope  with  a  positive  temperature  gradient  introduces  a 
radical  change  into  the  nature  of  the  tide.  Insofar  as  it  is’ still 
legitimate  to  refer  the  tides  in  such  an  atmosphere  to  those  in  an  ocean 
of  equivalent  depth  H,  it  can  be  said  that  H  becomes  a  function  of 
the  period,  which  increases  indefinitely  with  the  period.  The  bearing 
of  these  results  on  the  resonance  theory  of  atmospheric  tides  is  discussed. 


INTRODUCTION 


The  fact  that  the  solar  semidiurnal  tide  in  the  atmosphere  is  about 
100  times  larger  than  the  equilibrium  value  led  Kelvin  to  the  conclusion 
that  the  atmosphere  possesses  a  free  period  of  tidal  oscillation  of 
12  solar  hours.  The  resonance  period  of  the  atmosphere  must  be  within 
4  minutes  of  half  a  solar  day  in  order  to  account  for  the  observed  ampli¬ 
fication  and  for  the  fact  that  the  moon,  whose  tidal  force  is  more  than 
twice  as  large  as  that  of  the  sun,  excites  a  barometric  oscillation  having 
an  amplitude  of  only  one-sixteenth  that  of  the  solar  wave.  If  the  vertical 
temperature  distribution  in  the  atmosphere,  assumed  horizontally  strati¬ 
fied,  is  known,  it  is  possible  to  compute  its  free  period  of  tidal  oscil¬ 
lation  following  a  method  due  to  the  work  of  G.  I.  Taylor  (reference  l). 
According  to  this  method  one  first  determines  the  speed  of  propagation  of 
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long  waves  .V  in  a  flat  atmosphere  having  the  same  vertical  temperature 
distribution.  The  resonance  period  or  periods  of  tidal  oscillations  of 
the  atmosphere  are  then  identical  with  those  of  an  ocean  of  depth  H 
enveloping  the  earth  and  are  determined  from  the  relation 

V  =  (gH)^/^  (1) 

The  dependence  of  the  free  period  of  the  ocean  on  its  depth  H  is 
known  from  classical  tidal  theory.  For  the  solar  semidiurnal  oscilla¬ 
tion  E  =  7. 87  kilometers  ±  3  percent,  and  therefore  V  must  be 
0.278  kilometer  per  second  or  0.8l7co,  where  Cq  denotes  the  velocity 
of  sound  in  the  air  next  to  the  ground,  that  is,  3^  meters  per  second 
for  an  assumed  surface  temperature  of  288  K.  On  assuming  a  factor 
exp  i(fft  -  kx)  the  divergence  X  of  the  wave  motion  is  in  the  limit  of 
long  waves  found  to  satisfy  the  equation 


2  d^X 
c  - 

dz2 


+ 


^  +  g(7  -  1) 

/  dz  H 

dz 

=  0 


c^  =  7RT 
H  =  a2/gk2 
V  =  q/k 


(2) 


subject  to  the  boundary  conditions  that  at  the  ground  the  vertical  com- 
•  ponent  of  the  velocity  w  vanishes: 

W  =  ^^x/g)  +  nj^rx  +  (c2/g)(dx/dz^  =  0  (3) 

and  that  the  energy  of  wave  motion  per  unit  column  of  the  atmosphere  is 
not  infinite.  If  p(z)  denotes  the  density  of  the  air  as  a  function 
of  height  z  then  of  the  two  independent  solutions  of  equation  (2)  the 
one  is  accepted  which  gives  a  finite  value  for  the  energy  integral  U 


U  =  (1/2)  /  p(z) 
"^0 


+  V' 


w2) 


dz 


where  u,  v,  and  w  denote  the  components  of  velocity  due  to  the 

wave  motion.  » 
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When  the  temperature  at  great  heights  is  assumed  to  he  sufficiently 
low,  or  decreasing^  it  is  actually  found  (reference  2)  that  one  of  the 
solutions  of  equation  (2)  gives  a  finite  value  for  U,  while  in  the 
other  the  distribution  of  (u^  +  v^  +  w^)  with  height  is  such  that  U 
becomes  infinite  as  the  integral  is  extended  to  infinity.  When,  however, 
the  temperature  at  great  heights  is  assumed  to  increase  linearly  without 
limit,  it  is  found  that  for  both  solutions  of  equation  (2)  U  becomes 
infinite,  so  that  the  question  of  the  choice  between  the  two  solutions 
of  equation  (2)  is  left  open.  This  is  shown  in  a  following  section. 

This  work  was  conducted  under  the  sponsorship  and  with  the  financial 
assistance  of  the  National  Advisory  Committee  for  Aeronautics. 


SYMBOLS 


A  constant  in  equation  (63) 

A,  J,  B,  Y  coefficients  in  equations  (8)  and  (9) 

b  coefficient  in  equation  (52) 

C  =  sin  Jtb/sin  Jt(5  +  n) 

c  velocity  of  sound  at  any  level 

c^  velocity  of  sound  in  air  next  to  ground 


D  =  (2/7)  X  Ql  + 

E  =  (8]t/g)(l  +  n)cQQl  + 


F 

F3,  F4,  N 
F(0) 


H 


confluent  hypergeometric  function 
defined  by  equations  (33) 
defined  by  equations  (3^a) 

Airy  function  (reference  3) 
acceleration  of  gravity 
depth  in  a  uniform  ocean 


I 


Bessel  function  with  imaginary  argument 


k 
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J  Bessel  function  with  real  argument 

K  =  r(l  +  6)  sin  Jt(6  +  n)/2’^'*'^  sin  Jtn 
k  wave  number 

M  coefficient  in  equation  (47) 

m  =  [7.+  (7  -  1)(1  +  n)]^7H 

n  =  (g/Rp)  -  1 


P 

Q(u) 

R 

T 

t 

U 


U,  V,  w 
V 

W(x) 

z 

X 


y  =  2^[Ex  =  2^^" 
^0 


P 

r 


period  of  wave 

Brunt *s  period  (reference  4) 

function  defined  by  equation  (54) 

gas  constant 

temperature 

time;  in  appendix  the  plane  t 
energy  integral 

components  of  velocity  due  to  wave  motion 
speed  of  propagation  of  long  waves;  phase  velocity 
Whittaker* s  function  (reference  5) 
defined  by  equations  (19)  and  (20) 

distances  in  x-  and  z-directions 
roots  of  equation  (35) 
value  of  X  at  ground 

level  above  which  temper atiore  increases  linearly 
with  height 

constant  in  equation  (5) 
function  in  equation  (I9) 
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7  .  ratio  of  specific  heats 

6  =  -1  -  (n/2)  +  [t(1  +  n)/2:;^  +  (a)/27T) 

6  function  in  equation  (46) 

e  =  (1/6)8"^/^  1 48  -  xp/^ 

P  density  of  air 

pQ  undisturbed  density 

cr  frequency  of  wave 

T  =  a^/gk  =  kH 

roots  T  of  equation  (16) 

=  2/mz  -  (nJt/2)  -  (3V^) 
cp2  =  2^mz  +  (nJr/2)  +  {it/k) 

X  divergence  of  wave  motion 

(D  =  7  +  (7  -  l)(l  +  n) 

APPARENT  FAILURE  OF  ENERGY- INTEGRAL  CRITERION  FOR  AN  ATMOSPHERE 
IN  WHICH  TEMPERATURE  INCREASES  LINEARLY  WITH  HEIGHT 
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Equation  (2)  then  takes  on  the  form 

^d2x/dz2^  -  (n/z)(dx/dz) 

m  =  1/  +  (7  -  1)(1 

whose  solution  is 

X  =  Az(^‘'’’^)/^J]^^j^(2'JEz)  +  Bz(^'*‘^^^/%j^^^(2\Imz)  (8) 

when  n  is  an  integer,  or 

X  =  +  Bz(1+^)/2j_^_^(2\^)  (9) 


+  (m/z)x  =  0 
+  n^y^XH 


(7) 


when  n 


is  not  an  integer. 


In  evaluating  the  energy  integral  U  in  equation  (4)  for  the  pxirpose 
of  making  the  choice  between  the  solutions  in  equation  (8)  or  in  equa¬ 
tion  (9),  one  has,  for  a  wave  progressing  in  the  x-direction. 


u  = 

V  = 


(10) 


Since  c^X  RJ  zX  while,  for  large  values  of  z,  dx/dz  x  X  it  follows 
that  at  great  heights  both  u  and  w  vary  like  zX.  Now  for  large 
values  of  z 


cos  cp^  (ll) 

9^  =  j^2\^  -  (nit/2)  -  (3V**-)J 

J  ^  ~  (jt2mz)-lA  cos  cp2  (12) 

Cpg  =  j^2\(mz  +  (nJt/2)  +  (’^A)^ 
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It  follows  that  4-  +  w^)  in  the  integrand  of  equation  (4)  varies 

like  zl/2  cos^^  or  zl/2 

cos^2  the  two  solutions  of  equa¬ 

tion  (9)^  respectively.  The  kinetic  energy  per  unit  volume  therefore 

increases  ultimately  with  height  like  zl/2  and  the  energy  integral  U 
does  not  converge  for  either  solution.  .  Not  only  does  the  energy-integral 
criterion  not  enable  one  to  make  the  choice  between  the  two  solutions  of 
equation  (2)  but  also  neither  solution  yields  a  finite  value  for  U. 


RESOLUTION  OF  DIFFICULT!  WITS  THE  ENERGY  INTEGRAL 


The  result  obtained  in  the  previous  section  that  for  an  atmosphere 
in  which  the  temperature  increases  linearly  with  height  the  energy  inte¬ 
gral  of  the  wave  motion  diverges  for  both  solutions  of  the  wave  equa¬ 
tion  is  hard  to  reconcile  with  the  following  physical  considerations. 

A  mediijm  such  as  the  atmosphere^  which  is  bounded  by  a  perfectly 
reflecting  surface  at  the  ground  but  is  exposed  to  space  on  the  other 
side,  would  be  expected  to  be  facilitated  in  sustaining  a  free  oscilla¬ 
tion  the  more  it  is  capable, of  diverting  wave  energy  from  great  heights 
toward  the  ground.  The  more  the  upper  atmosphere  behaves  like  a  reflector, 
the  greater  will  be  its  trapping  power  for  waves,  and,  with  it,  its 
ability  to  propagate  waves  to  great  ranges  -  a  characteristic  of  free 
oscillations.  Now  it  is  known  (references  6  to  8)  that  an  atmosphere  of 
constant  temperature  throughout  can  completely  trap  waves  either  in  the 
whole  spectrum  or  in  a  limited  frequency  range.  This  is  also  true  for 
an  atmosphere  of  the  Taylor  type  in  which  the  temperature  is  assumed  to 
decrease  at  a  constant  rate  in  the  troposphere  and  to  remain  constant  in 
the  stratosphere.  Now  when  the  temperature  is  constant,  the  rays  ema¬ 
nating  from  a  source  situated  in  the  atmosphere  are  rectilinear,  and  all 
but  the  horizontal  rays  are  lost  to  space.  In  the  Taylor  type  atmosphere 
the  rays  are  even  bent  upward  in  the  troposphere,  eventually  to  become 
rectilinear  in  the  stratosphere.  In  spite  of  these  nonconservative  prop¬ 
erties,  both  types  of  atmosphere  are  capable  of  sustaining  free  oscil¬ 
lations  in  certain  frequency  ranges.  In  the  case  of  an  atmosphere  with 
a  linearly  increasing  temperature  all  the  rays  are  bent  downward,  so 
that  its  trapping  powers  should  be  even  greater.  Yet,  the  divergence  of 
the  energy  integral  found  in  the  previous  section  indicates  a  tendency 
to  allow  the  waves  to  leak  out  to  higher  elevations. 

The  resolution  of  this  paradox  lies  in  the  fact  that  equation  (2) 
for  the  divergence  X  is  not  exact  but  is  an  approximation  for  the 
limiting  case  of  long  waves.  This  approximation  breaks  down  when  the 
temperature  increases  indefinitely  linearly  with  height.  The  exact 
equation  for  X  is  equation  (13) 


+  fi  doi  _  ZfiW  J 

- 1)  + 

1  1  dc2  g(7  -  1) 

1 

dz2  Vc2  dz  c2/  <iz 

H  ,2  dz  g,2 

J 

x=  0  (13) 
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with 


(1  -  H2k2)w  =  (c^x/g)  +jE  -7X  +  (c^/g)  (dx/dz^  (l^) 

replacing  equation  (3).  Now  in  the  limit  of  long  waves,  when  the  wave 
number  k  vanishes,  the  long-wave  approximation  is  made  by  dropping 
the  k^-term  in  equation  (13).  This  approximation  is  valid  in  an  atmos¬ 
phere  of  constant  or  decreasing  temperature.  In  the  case  of  an  atmos¬ 
phere  in  which  the  temperature  (and  with  it  c^)  increases  linearly  with 
height,  the  terms  in  brackets  in  equation  (13)  decrease  like  l/z  and 
eventually  become  smaller  than  the  k^-term.  In  fact  at  extremely  large 
heights  equation  (13)  reduces,  for  a  given  k,  to 


(d2x/dz2)  -  k2x.=  0  (15) 

indicating  that  the  appropriate  solution  for  X  ultimately  varies  like 
e“^^  and  not  according  to  the  power  law  that  would  follow  from  equa¬ 
tions  (9),  (11),  and  (12).  When  the  complete  equation  (13)  is  used  one 
finds  that  the  solution  having  the  factor  e"^^  yields  a  finite  value 
for  the  energy  integral  U  in  equation  (^),  whereas  the  other  solution, 
with  the  factor  e^,  must  be  discarded  because  it  gives  a  divergent 
value  for  U.  With  the  complete  wave  equation  (13)  there  is  therefore 
no  ambiguity  as  to  the  choice  between  the  two  independent  solutions,  nor 
is  there  any  difficulty  with  the  energy  integral. 

Using  the  notation  of  equations  (5)  and  putting 


X  =  2kz 

T  =  cT^/gk  =  kH 

5  =  _1  _  (n/2)  +  [t(1  +  n)/2^  +  (cd/27t) 
CD  =  7  +  (7  -  1)(1  +  n) 


(16) 


equation  (13)  reduces  to 


(d^x/dx^)  -  n(dx/dx)  +  +  1  +  (n/2)  -  (x/l^)^X  =  0  (l?) 
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Choosing  the  solution  of  equation  (17)  which  is  bounded  at  infinity, 

X  =  x^/%(x)  .  .  — ^  il  - 

(2S+24-n)  (1+n) 

2^2  (18) 

where  W(x)  is  Whittaker's  fimction  (reference  5,  p.  337).  The  expo¬ 
nential  behavior  of  X  for  large  values  of  x  assures  the  convergence 
of  integral  (4)  for  U,  whereas  in  the  second  solution  of  equation  (17) ^ 

which  behaves  like  x"^"^e^/^  at  infinity,  the  integral  diverges. 


LIMITING  FORM  OF  SOLUTION  FOR  LONG  WAVES 


j6(5  +  1  +  n)/:^  +  .  .  J- 


With  the  solution  of  the  exact  equation  (13)  given  by  equation  (I8) 
it  becomes  of  interest  to  investigate  the  form  it  assmes  in  the  limit 
of  long  wavelengths.  In  particular  the  relation  of  this  limiting  form 
of  equation  (I8)  to  solution  (9)  of  the  approximate  equation  (2)  is  to 
be  investigated,  with  a  view  of  throwing  light  on  the  question  of  the 
choice  of  the  proper  linear  combination  of  the  two  independent  solutions 
appearing  there.  For  this  purpose  the  following  relations  (reference  5, 
PP.  337^  338,  and  3^8)  are  used; 


where  F  denotes  the  confluent  hypergeometric  fxmction.  Thus,  for  the 
case  when  n  is  not  an  integer  there  is  obtained 

X(x)  =  [r(-l  -  n)/r(-5  -  n  -  ij]  xl+ne-^/%(-6,  2  +  n,  x)  + 

[r(l  +  n)/r(-5^e-x/2F(-5  -  n  -  1  -  n,  x) 


(21) 
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X(x)  =  |[r(2  +  5  +  n)  sin  «(&  +  nj]/(r(2  +  n)  sin  xl+ne-V%(-5,  2  + 

n,  x)  -  |r(l  +  n)  X  r(l  +  &)  sin  (n8)/j^e"^/2F(-5  -  1  -  n,  -n,  x) 

(22) 

where  in  the  last  transformation  the  following  relation  was  used 


r(x)  X  I*(l  -  x)  =  «/sin  rtx  (23) 

Passing  now  to  the  limit  of  long  waves, 

T  =  KB  «  1 
5— >a5/27T  =  m/2k  »  1 

With  8  very  large,  the  F-functions  in  equation  (22)  reduce  to  Bessel 
functions: 

F(-8,  2  +  n,  x) — »jl  -  8x/(2  +  n)  +  (8x)^/2(2  +  n)(3  +  n)  +  •  •  = 

r(2  +  n)(5x)-(^+l)/2j^+l(2\l6x)  (25) 

F(-5  -  1  -  n,  -n,  x) — >r(-n)(8x  )(n+l)/2j____l{2^)  (26, 

It  is  important  to  observe  that  these  limiting  forms  of  the  F*s  are 
valid  only  when  both  6  »  1  and  x  «  b,  for  otherwise  the  series 
for  F  depends  on  the  higher  powers  of  x,  when  it  is  no  longer  legiti¬ 
mate  to  replace  the  factor  S  -  k  by  5  in  the  coefficients.  Sub¬ 
stituting  now  equations  (25)  and  (26)  into  equation  (22)  and  performing 
some  further  reductions  with  the  aid  of  equation  (23)^  there  is  obtained 


x  =  K  exp|I(y2/88)^  jy’^+Vn+l(y)  +  Cy’^'^^J_n_i(y^ 


(27) 
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y  =  2^Ex  =  2^^z 

K  =  r(l  4*  6)  sin  Jt(6  +  n)/2^'^^  sin  nn 
C  =  sin  Tr5/sin  jr(&  +  n) 


(28) 

(29) 


Except  for  the  factor  exp  which  is  slowly  variable  in  the 

limit  of  large  values  of  5,  equation  (27)  is  identical  with  equation  (9)* 
The  ratio  of  constants  B/a  in  the  latter  is  now^  however,  no  longer 
arbitrary  but  is  a  definite  function  C  of  5,  as  given  in  equation  (29). 


The  goal  of  determining  the  appropriate  solution  of  the  wave  equation 
for  an  atmosphere  with  a  linear  increase  of  temperature  has  thus  been 
achieved,  but  on  closer  study  the  physical  implications  of  the  results 
are  rather  disturbing.  In  order  to  ascertain  the  physical  significance 
of  the  parameter  6  the  phase  velocity  and  the  period  are  determined. 

The  phase  velocity  V  is  given  by 


=  (^/kcQ  =  {2I7)  X  |Tl  +  n)(q^^^(l/y)  =  D/y  (30) 

while  the  period  of  the  wave  P,  measured  in  seconds,  is  given  by 

P  =  2n/a  =  (8rt/g)(l  +  n)co[(l  +  n)c^  =  E5/y  (31) 

For  n  =  U.5,  then  D  =  6.36  and  E  =  IO8O,  while  for  n  =  8.5,  which 
is  appropriate  for  the  positive  lapse  rate  prevailing  in  the  NACA  model 
atmosphere  between  83  and  120  kilometers  (reference  9),  B  =  10.0  and 
E  =  1180.  For  v/cQ  of  the  order  of  unity,  the  period  P  changes  by 
about  3  minutes  as  5  changes  by  unity.  But,  since  C(6)  in  equation  (29) 
is  a  periodic  function  of  6  of  period  unity  and  takes  on  all  values 
between  -00  and  »  as  6  changes  by  1,  it  follows  that  in  the  limit  of 
long  periods  X  does  nob  approach  a  fixed  form  but  changes  its  character 
radically  as  the  period  is  changed  by  only  a  few  minutes.  This  increased 
sensitivity  of  X  to  the  period  for  an  atmosphere  with  positive  lapse 
rate  is  to  be  contrasted  with  the  decreasing  sensitivity  to  change's  in 
period  for  atmospheres  with  zero  or  negative  lapse  rates. 
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FREE  OSCILLATIONS  OF  AN  ATMOSPHERE  IN  WHICH  TEMPERATURE 
INCREASES  LINEARLY  WITH  HEIGHT 


In  order  to  gain  more  insight  into  the  acoustic  propagating  prop¬ 
erties  of  an  atmosphere  with  a  constant  positive  tenperature  gradient, 
the  free  oscillations  or  such  an  atmosphere  when  it  extends  from  the 
ground  up  will  he  investigated.  Therefore  the  variation  of  the  phase 
velocity  V  with  the  period  P  for  the  first  and  for  the  higher  normal 
modes  will  he  studied.  Such  a  variation  is  shown,  for  example,  hy 
curve  I  in  figure  2,  for  the  first  mode  of  an  atmosphere  with  constant 
negative  temperature  gradient  (references  8  and  10 ) .  For  a  Taylor  type 
atmosphere  having  the  same  negative  temperature  gradient  in  the  tropo¬ 
sphere,  the  corresponding  curve  starts  with  a  value  of  0.86  at 
P  =  2  minutes  and  reaches  its  limiting  value  of  about  0.92  at 
P  =  4  minutes  (references  8  and  10). 

The  divergence  X  from  equation  (22)  can  he  written  in  the  form 

X  =  xl+ne-x/2j.3^  _  Ge'^/^Fp  (32) 


with 


G  =  [r(l  +  n)r(2+n)r(l  +  6)  sin  JtnXsin  Jt0y^[rtr(2  +  S  +  n) X  sin  )t(8 +n)] 


Ft  =  F(-5,  2  +  n,  x) 


Fg  =  F(-5  -  1  -  n,  -n,  x) 


The  possible  free  oscillations  of  the  atmosphere  are  determined  from  the 
condition  of  the  vanishing  of  the  vertical  component  of  velocity  w  at 
the  ground,  which,  according  to  equation  (l4),  is  given  hy 

(1  +  n)(l  -  t2)(w/7HX)  =  -1  -  n  +  (x/2t)  +  (x/x)(dX/dx)  =  0  (3^) 


With 


Fj  =  F(l  -  6^  3  +  n,  x) 

F4  =  F(-8  n,  1  -  n,  x) 

D  =  F]^  -  Gx"^"%2 

N  =  [(1  +  n)/:^Fi  -  [8/(2  +  n[[F3  -  g[(8  +  1  +  n)/^x-l-“F4 
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equation  (3^)  leads  to  the  condition 

(1  +  n)(l  -  t^VwAhx)  =  -(l  +  n)  +  (x/2)  [(i/t)  -  l”]  +  xN/d  =  0 

^  ^  I  ^35) 

X  =  Xq  J 

where  xq  denotes  the  value  of  x  at  the  groimd. 

The  problem  of  determining  the  dependence  of  the  phase  velocity  V 
on  the  period  P  involves  the  determination  of  the  dependence  of  the 
roots  X  of  equation  (35)  on  the  parameter  &.  For  a  given  5  one 
first  determines  the  two  roots  t  of  equation  (l6):  . 

=  ^7(8  +  1  +  (n/2^  +  «  7^|5  +  1  +  (n/2^  -  (1  +  n)a>^  ^  ^ 

T2  =  frfs  +  1  +  (n/2^  -  <;  +  1  +  (n/2y  -  (l  +  n)a3^  ^^1  +  n) 

^  (37) 

For  each  t  one  finds  the  roots  x  of  equation  (35)  end  determines  V 
and  P  from 


V/cq  =  [2(1  +  n)T/7x] 


(38) 

n)/7T3d 

(39) 

The  free  oscillations  of  an  atmosphere  have  heen  determined  with 
n  =  ^1.5^  corresponding  to  a  positive  temperature  gradient  of  about 
6.4°  per  kilometer.  Figure  1  shows  the  variation  of  the  roots  x  with 
6  for  the  two  branches  and  T2.  It  is  shown  in  appendix  A  that 

for  the  first  mode  in  the  which  is  shown  by  curve  A  in  fig¬ 

ure  1,  the  lower  end  of  the  curve  is  given  by 

X  =  1.88(1  -  8)^/^-5  (to) 

for  sufficiently  small  values  of  x.  With  T2  =  follows  from 

equations  (38)  and  (39)  that  along  this  end  of  curve  A,  p— >00^ 

V/cq“->w,  and  v/cq->(i/10)P^  when  P  is  measured  in  minutes.  Hence 
in  an  atmosphere  with  constant  positive  temperature  gradient,  the  phase 
velocity  does  not  approach  a  limit  with  increasing  period  but  increases 
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linearly  with  the  period.  The  same  result  applies  to  the  higher  modes. 
Figure  2  shows  the  ftinction  V(P)  which  results  from  the  x(6)-curves  in 
figure  1.  It  is  shown  in  appendix  B  that  in  the  Tg-hranch  no  waves  can 
he  propagated  for  periods  less  than  4.5  minutes,  which  is  Brunt's  period 
(reference  4)  determined  from 

=  g^(7  -  1)  +  g(dc^/dz^  (4l) 

In  all  the  modes  v/cq  grows  indefinitely  with  the  period. 

It  is  of  interest  to  see  the  type  of  dispersion  one  gets  hy  using 
the  "long  wave"  solution  (27).  Dropping  the  exponential  term. 


X  =  y’"'-ljn+i(y)  +  Cy’"+lj_n_i(y) 

(42) 

dx/dy  =  y’^'^^JnCy)  -  Cy’^'’’^J.n(y) 

(43) 

while  the  seculsir  equation  (35)  reduces  to 

-(1  +  n)  +  (7y^/4n)  +  (y/2)  |jn(y)  -  CJ_n(y)J^|jn+l(y)  +  CJ_n-l(y^  =  0 

(44) 


In  this  case  one  can  solve  directly  for  C,  since  all  the  other  terms 
do  not  depend  explicitly  on  6; 


C  = 


-  jTl  +  n)(2/y)  -  (7y/2aj^Jn+i  +  n)(2/y) 

(7y/2aDQj_n_i(y)^  (^5) 


With  C(5)  given  in  equation  (29),  equation  (45)  allows  one  to  obtain  5 
as  a  function  of  y,  from  which  one  determines  the  function  V(P)  hy 
using  equations  (30)  and  (31).  The  results  are  shown  in  figure  3- 


It  is  seen  that  the  general  shape  of  the  V(P) -curves  is  similar  to 
that  of  the  exact  ciirves  given  in  figure  2.  The  deviations  are  due  to 
the  fact  that,  in  the  upper  right-hand  corner,  8  is  of  the  order  of 
unity,  for  which  the  long- wave  approximation  cannot  apply.  The  curves 
in  figure  3  become  also  inaccurate  as  the  cut-off  period  is  approached. 
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because  there  y  is  of  the  order  of  5,  and  the  Bessel  function  approxi¬ 
mation  in  relations  (25)  and  (26)  breaks  down.  The  analysis  for  this 
limiting  case  is  given  in  appendix  B. 


APPLICATION  TO  THEORY  OF  ATMOSPHERIC  TIDES 


This  investigation  has  shown  that  in  an  atmosphere  where  the  temper¬ 
ature  increases  indefinitely  with  height  at  a  constant  rate,  the  acoustic 
propagation  properties  are  radically  different  from  those  that  have 
hitherto  been  encountered  in  the  theory  of  atmospheric  tides.  The 
principal  new  feature  is  that  the  phase  velocity  V  no  longer  reaches 
a  limit  for  waves  of  long  period.  For  atmospheres  where  V  reaches  a 
limit,  the  theory  of  atmospheric  tides  is  simplified  by  an  application 
of  Taylor's  theorem  (reference  1),  according  to  which  the  tides  in  the^ 
atmosphere  are  equivalent  to  those  in  an  ocean  of  uniform  depth  H  =  V  /g. 
In  our  case,  V  grows  indefinitely  with  the  period,  and,  insofar  as  it 
is  still  legitimate  to  refer  to  an  equivalent  ocean  of  depth  H,  one 
may  say  that  H  depends  on  the  period,  ultimately  increasing  as  the 
square  of  the  period.  What  happens  is  that,  as  the  period  is  increased, 
the  distribution  of  energy  of  wave  motion  with  height  no  longer  approaches 
a  fixed  form,  but  the  energy  continues  to  spread  to  higher  elevations  as 
the  period  is  increased.  A  discussion  will  not  be  given  here  of  the 
nat\ire  of  the  tides  £n  an  atmosphere  with  constant  positive  temperature 
gradient  extending  from  the  ground  up.  In  such  an  atmosphere  the  tides 
are  probably  different  from  those  in  an  ocean  of  fixed  depth.  An  investi¬ 
gation  of  this  problem  would  be  in  order  if  the  evidence  would  continue 
to  accumulate  to  the  effect  that  the  temperature  in  the  ionosphere 
increases  continuously  with  height. 

Pending  such  an  investigation,  suffice  it  to  remark  that  there  is 
no  difficulty  in  treating  the  case  where  a  layer  with  positive  tempera¬ 
ture  gradient  is  topped  by  an  outer  envelope  of  constant  or  decreasing 
temperature.  In  that  case  the  proper  solution  for  the  layer  is  given 
by  equation  (9),  with  B/a  arbitrary.  It  would  also  follow  that  in  the 
absence  of  proof  as  to  the  temperature  of  the  outer  envelope  of  the  atmos¬ 
phere  one  should  not  in  theoretical  studies  allow  an  assumed  positive 
gradient  in  an  interior  layer  to  extend  through  the  outer  envelope.  The 
consequences  of  such  an  extrapolation  of  the  temperature  curve,  if  rigor¬ 
ously  determined,  ar*e,  in  the  light  of  the  results  of  this  investigation, 
likely  to  affect  seriously  the  natvire  of  the  atmospheric  tides.  Any 
preliminary  results  obtained  for  an  atmosphere  with  a  positive  tempera¬ 
ture  gradient  extending  through  the  envelope  should  therefore  be  checked 
by  comparing  them  with  those  for  a  modified  atmosphere  in  which,  say, 
the  temperature  in  the  outer  envelope  is  assumed  to  be  constant  or 
decreasing  with  height.  It  would  also  be  helpful  to  compute  in  each 
case  the  distribution  with  height  of  the  energy  density  of  the  wave 
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motion  (p/2)^^  +  v2  +  w2)  in  order  to  assure  that  any  segment  of  the 
temperature  curve  which  is  as  yet  not  supported  by  observations  does  not 
effect  a  marked  redistribution  in  the  wave  energy. 


Institute  for  Advanced  Study 

Princeton,  N.  J-,  September  26,  19^9 
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APPEMDIX  A 

V 

DEPENDENCE  OF  PHASE  VELOCITY  V  ON  PERIOD  P  WHEN 
PERIOD  IS  INCREASED  INDEFINITELY 


Only  the  r^-hranch  of  the  first  mode  will  be  treated  here  and  in  the 

particular  case  when  n  =  4.5*  The  higher-order  modes^  as  well  as  those 
of  the  T2_ -branch,  can  be  dealt  with  similarly.  The  problem  is  to  determine 
the  limiting  form  of  the  root  x(5)  of  equation  (35)  as  5—^1  -  €.  Then 


■^2  =  Vll  ] 

(i^6) 

G-»  -  6r(5.5)/6.5  J 

Also  for  a  vanishing  x  all  the  F’s  reduce  to  unity. 

Now  assume  that 

G— 

(47) 

as  X— ^0,  and  seek  to  determine  the  coefficient  M  from  equation  (35) 
On  retaining  only  leading  terms,  equation  (35)  becomes 

-5.5  +  0.875X  +  [5.5  -  (x/6.5^(1  +  Mx)  =0 

■  (48) 

On  annulling  the  coefficient  of  x  there  is  obtained 

M  =  -0.721  =  -(1  -  5)|r(5.5)/6.^x“^-5 

(49) 

X  =  1.88(1  -  6)V6-5 

(50) 
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APPENDIX  B 

DETERMINATION  OF  CUT-OFF  PERIOD 


If  equation  (^5)  is  evaluated  for  large  values  of  y  By  using  the 
asymptotic  forms  of  the  Bessel  functions  it  is  found  that  all  modes 
coalesce  at  a  limiting  period  (5.7  min  for  n  =  ^+.5)^  for  which  the 
phase  velocity  vanishes.  This  result  can,  however,  not  he  relied  upon 
since  when  y  becomes  of  the  order  of  &  the  Bessel  function  approxi¬ 
mations  in  relations  (25)  and  (26)  are  no  longer  valid.  In  order  to 
determine  the  cut-off  period  the  asymptotic  form  of  W  in  equation  (18) 
has  to  be  obtained  when  both  6  and  x  become  large  and  are  of  the 
same  order  of  magnitude.  For  this  purpose  Whittaker's  (reference  5, 
p.  339)  countour- integral  representation  for  W  is  used,  which  can 
be  put  in  the  form 

X  =  xV2w(x) 

=  (l/2Kl)ef’^r(l  +  6)e-V2xl+n^t-S-l(l  +  dt  (51) 

the  path  starting  at  «»  and  returning  to  «>  after  encircling  the  origin 
in  the  positive  direction,  but  excluding  the  point  t  =  -1.  Now  put 


X  =  45  -  b5 
t  =  u  -  (1/2) 
then  near  t  =  -(l/2)  there  is  the  expansion 


(52) 


t-&-l(i  +  t)S+l+ne-xt  ^ 


^  ejqjj^i’fS  +  (x/2)  +  Q(u^  (53) 


Q(u)  =  (6b)u  +  (165/3)u3  +  (61(-6/5)u5  +  .  .  . 


O 


(l/2J(i)r(l  +  6)x^‘'‘’^^ 


u 


exp  Q(u)  du 


(54) 

(55) 


where  the  path  in  the  u- plane  is  similar  to  the  former  path  in  the  .t-plane 
For  large  values  of  6  and  for  values  of  b  of  the  order  of  S“  /j. 
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the  principal  contribution  to  the  integral  in  equation  (55)  comes  from 
the  immediate  vicinity  of  the  origin.  This  can  be  made  evident  by  putting 


(i66/3)u3  =  v3 
A  =  (3/i6)^/1)5^/3 


(56) 


when  the  u^-term  in  equation  (5^)  becomes  of  the  order  of  while 
the  expression  in  brackets  in  the  integrand  of  equation  (55)  reduces  to 
-2”^  to  within  terms  of  the  order  of  6-1/3. 

Hence, 


exp  ^v3  -f  Av^  dv  (57) 

Here  the  path  starts  at  oo  and  terminates  at  ooe'  -iTt/3  after 

encircling  the  origin  in  the  positive  direction.  The  integral  can  be 
evaluated  in  terms  of  the  Airy  f-unction  (reference  3>  P«  l89)>  yielding 


X  Ri -(l/2«i)2-’^xl+’^r(l  +  6) 


X  R*  2-'^xl+'^r(l  +  6)(2886)-1/3f(0) 

0  =  (1/6)5-1/2|46  -  x|3/2 


(58) 


F(0)  =  0l/3|j_^^^(0)  + 


<  46 


(59) 


F(e)  =  0i/3ji_^^^(0)  -  i^/3(e^ 


> 


X  >  45 


(60) 
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Equation  (6o)  shows  that  heyond  !)•&,  X(x)  is  no  longer  oscillatory. 
The  largest  zeros  of  X  are  determined  from  the  zeros  of  F(0)  (refer¬ 
ence  3,  p.  751)  in  equation  (59): 


0^  =  2.383 
=  5.610 
0^  =  8.657 


(61) 


On  substituting  expression  (59)  into  the  secular  equation  (35)  one  finds 
that  in  the  limit  of  large  values  of  6  the  roots  of  equation  (35) 
coincide  with  the  roots  of  X-  Hence  the  largest  roots  of  equation  (35) 
are  given  by 


Xi  =  46  -  ^3650^^) 


1/3 


(62) 


where  the  0i's  are  given  in  equation  (6I).  From  equation  (62)  it 
follows  that  near  the  cut-off  period 


V/cq— ^(a/&)  |l  +  (l/8)(60i/&)2/^ 
P  =  PbIT  +  (l/8)(60i/8)2/^ 


(63) 


where  A  is  a  constant  and  Pg  denotes  Brunt's  period  (reference  4). 
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Figure  1.-  Dependence  of  roots  x  of  equation  (35)  on  8  for  an  atmosphere 
with  constant  positive  temperature  gradient,  n  =  4.5.  Solid  curves  for 

dashed  curves  for  t  .  Encircled  numbers  refer  to  order  of  mode. 
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an  atmosphere  with  constant  positive  temperature  gradient,  n  =  4.5.  The  -modes  start  at  Brunt’s 

period  of  4.5  minutes.  Curve  I  is  for  an  atmosphere  with  constant  negative  temperature  gradient  of 
seven-elevenths  of  the  adiabatic.  Encircled  numbers  refer  to  order  of  mode. 
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refer  to  order  of  mode. 


